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• Premise of the study: Microsatellite loci were developed in Sebaea aurea (Gentianaceae) to investigate the functional role of 
diplostigmaty (i.e., the presence of additional stigmas along the style). 

• Methods and Results: One hundred seventy-four and 180 microsatellite loci were isolated through 454 shotgun sequencing of 
genomic and microsatellite-enriched DNA libraries, respectively. Sixteen polymorphic microsatellite loci were characterized, 
and 12 of them were selected to genotype individuals from two populations. Microsatellite amplification was conducted in two 
multiplex groups, each containing six microsatellite loci. Cross-species amplification was tested in seven other species of 
Sebaea. The 12 novel microsatellite loci amplified only in the two most closely related species to S. aurea (i.e., S. ambigua and 
S. minutiflora) and were also polymorphic in these two species. 

• Conclusions: These results demonstrate the usefulness of this set of newly developed microsatellite loci to investigate the mat- 
ing system and population genetic structure in S. aurea and related species. 
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Sebaea aurea (L. f.) Roem. & Schult. is a diploid (2n = 28; 
Kissling et al., 2008) annual herb found in the Western Cape 
(South Africa). Sebaea aurea is a diplostigmatic species, with 
styles showing two stigmatic areas, one each at the apex and 
base of the styles (see Kissling et al., 2009a). This character is 
suspected to provide reproductive assurance and reduce seed 
discounting (i.e., the formation of self-fertilized seeds from 
ovules that, if they had not been self-fertilized, would have been 
cross-fertilized) and, therefore, might provide some evolutionary 
advantages (Kissling et al., 2009b). 

To date, no polymorphic microsatellite markers have been 
developed in this species or in other species for which micro- 
satellites could be potentially transferred to S. aurea. Therefore, 
the purpose of this study was to isolate and characterize micro- 
satellite loci in S. aurea to conduct reproductive biology experi- 
ments and to estimate the selfing rate in a progeny array, which 
will contribute to a better understanding of the functional role 
of diplostigmaty. Furthermore, this set of microsatellite loci 
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will contribute to the investigation of the distribution of genetic 
diversity across the species range. 

METHODS AND RESULTS 

Microsatellite isolation was carried out through two independent 454 
pyrosequencing strategies. First, DNA was extracted from 20 individuals of 
S. aurea from the Helderberg Nature Reserve population (Appendix 1) using 
the DNeasy Plant Mini Kit (QIAGEN, Hombrechtikon, Switzerland) following 
the manufacturer instructions and then pooled and concentrated. Five hundred 
nanograms of total DNA was used to construct 454 genomic libraries by the 
sequencing service from the University of Valencia (Servicio Central de 
Soporte a la Investigation Experimental [SCSIE], Valencia, Spain) and shot- 
gun sequenced on a GS Junior 454 sequencer (454 Life Sciences, a Roche 
Company, Barcelona, Spain). We obtained 120,157 reads with an average read 
length of 468.43 bp that were trimmed of adapter and low-quality regions and 
assembled into contigs using GS De Novo Assembler implemented in Newbler 
2.5pl (454 Life Sciences, a Roche Company, Madrid, Spain). Generated con- 
tigs and unique reads not assigned to contigs were subjected to BLAST analy- 
sis, and those matching organellar (chloroplast or mitochondria) sequences 
were discarded. 

We screened all 64,422 unique reads and nonorganellar contigs with iQDD 
version 1.3.0.0 software (Meglecz et al., 2010). We set the script to identify all 
possible di-, tri-, tetra-, penta-, and hexanucleotide repeats with a minimum of 
five repeat units, as well as compound repeats, and to directly design primers 
using Primer3 (Rozen and Skaletsky, 2000). After discarding reads with too 
short flanking sequences, primers were successfully designed for 174 reads. Of 
these 73, 90, and 11 corresponded to di-, tri-, and tetranucleotide microsatellites, 
respectively. 

Additionally, size- selected fragments from genomic DNA were enriched for 
microsatellites by Ecogenics GmbH (Zurich- Schlieren, Switzerland) using 
streptavidin-coated magnetic beads and biotin-labeled (CT) 13 and (GT) 13 repeat 
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probes. The simple sequence repeat (SSR)-enriched library was sequenced on a 
Roche 454 platform using the GS FLX Titanium reagents at Ecogenics GmbH. The 
34,992 reads sequenced had an average length of 427 bp. Of these, 996 contained a 
microsatellite insert with a tetra- or a trinucleotide motif of at least six repeat units 
or a dinucleotide motif of at least 10 repeat units. Suitable primer design using 
Primer3 (Rozen and Skaletsky, 2000) was possible in 180 reads. Thirty-nine and 36 
loci from the first and second sequencing strategies, respectively, were tested for 
functionality and polymorphism using at least seven individuals. 

Amplifications were carried out in 10-jiL reactions containing 1-3 ng of 
template DNA, lx GoTaq Flexi Buffer (Promega Corporation, NeucMtel, 
Switzerland), 1.5 mM MgCl 2 , 0.2 mM dNTPs, 0.04 uM forward primer with a 
5' Ml 3 tail, 0.16 uM 6-FAM-labeled Ml 3 primer, 0.16 jiM reverse primer, and 
0.5 U HotStarTaq DNA Polymerase (QIAGEN). PCRs were performed in a 
Dyad Cycler (Bio-Rad GmbH, Cressier, Neuchatel, Switzerland), and the PCR 
cycling conditions consisted of an initial activation step of 5 min at 95 °C; fol- 
lowed by 30 cycles each of 45 s at 95°C, 60 s at 55°C, and 60 s at 72°C; eight 
cycles each of 45 s at 95°C, 60 s at 53°C, and 60 s at 72°C; and a final extension 
step of 30 min at 72°C. PCR products were run on an ABI3730xl Genetic Ana- 
lyzer (Applied Biosystems, Rotkreuz, Switzerland) using GeneScan 500 LIZ as 
an internal size standard (Applied Biosystems). Amplified fragment lengths 
were assigned to allele sizes with GeneMapper (Applied Biosystems). This 
screening detected 16 polymorphic loci that consistently amplified (Table 1); 
however, only 12 with adequate sizes for multiplex amplification were selected 
to genotype the whole set of samples. 

The 12 selected microsatellite loci were arranged into two multiplex reac- 
tions (I and II; Table 1). Amplifications were carried out in a 10-jiL reaction 
using Type-it Microsatellite PCR kits (QIAGEN). The PCR cocktail included 
5.0 jllL 2x Type-it Master Mix, 1.0 jiL primer mix (2 jiM each forward fluo- 
rescent-labeled and unlabeled reverse primer; Table 1), 3.0 uL RNase-free 
water, and 1 JiL template DNA (10-50 ng). The PCR program followed the 
manufacturer's conditions: an initial activation step at 95°C for 5 min; 28 cycles 
at 95°C for 30 s, 57°C for 90 s, and 72°C for 30 s; followed by a final extension 
step at 60°C for 30 min. PCR products were run and sized as in the aforemen- 
tioned pilot study. 

Table 1. Characteristics of 16 microsatellite loci developed in Sebaea aurea. 



Genotypic data were obtained from three populations of S. aurea (Helderberg 
Nature Reserve, N = 38; Somerset West, N = 25; Paradyskloof, N = 30; Table 2, 
Appendix 1) for the 12 microsatellite loci (Tables 1 and 2). Number of alleles 
(A), observed heterozygosity (H 0 ), and unbiased expected heterozygosity (H e ) 
(Nei, 1978) were calculated with GENETIX version 4.05 (Belkhir et al., 2004). 
Linkage disequilibrium between pairs of microsatellite loci and between pairs 
of microsatellite loci and populations, using 10,000 permutations, were calcu- 
lated with GENEPOP version 4.1.4 software (Rousset, 2008). This same soft- 
ware was used to estimate the frequency and 95% confidence intervals of 
estimated null allele frequencies of each locus in each population (Dempster 
et al., 1977). Of the 186 available pairwise comparisons between loci and pop- 
ulations, only seven showed significant linkage disequilibrium (P < 0.05), 
whereas of the 66 available pairwise comparisons between loci only two 
showed significant linkage disequilibrium (P < 0.05). Of these, only the pair 
Sebal 1-Sebal3 was consistent between the Helderberg and Paradyskloof popu- 
lations and the global estimate of linkage disequilibrium; however, none of the 
pairwise comparisons were significant after Bonferroni correction. 

We detected a total of 164 different SSR alleles for the 12 polymorphic 
microsatellite loci in the three analyzed populations of S. aurea. The number of 
alleles ranged from a minimum of three alleles for locus Sebal 9849 to a maxi- 
mum of 29 alleles for locus Sebal 3, and the mean number of alleles per locus 
was 13.67 ± 7.13. One locus (Sebal9849) was monomorphic in the Somerset 
West population (Table 2). H 0 ranged from 0.0 (locus Seba28206) to 0.868 
(locus Sebal 3) in the Helderberg population, and unbiased H e in the polymor- 
phic loci ranged from 0.050 (locus Sebal9849) to 0.945 (locus Seba28206) 
(Table 2). Almost half of the loci in each population showed significant devia- 
tion from Hardy-Weinberg equilibrium toward heterozygote deficiency, 
whereas the remaining loci in each population showed either nonsignificant 
heterozygote deficiency or excess. Global population inbreeding estimates gave 
significant heterozygote deficiency in all three populations. A likely expla- 
nation for this result may be the presence of null alleles. In fact, 17 out of the 
36 loci per population comparison gave estimated frequencies of null alleles 
higher than 0.05. Nonetheless, local population substructure or high selling 
rates within populations could also result in high inbreeding coefficient (F IS ) 



Locus 




Primer sequences (5'-3') a 


Repeat motif 


Allele size (bp) 


r a (°o 


Multiplex group 


GenBank accession no. 


Seba04 


F 
R 


VIC-ATTCACATCGCTTACAGCCC 
AATTTAAGAACGTCGCCGC 


(AG) 10 


130 


56 


I 


KF218835 


SebalO 


F 
R 


PET-CCTTTATGTGCAACGGGAAG 
AATGATCCATGCATTCTGCC 


(CT) 8 


135 


56 


I 


KF2 18836 


Sebal 1 


F 
R 


6FAM-TGCCTGTTCGACTTATGCAA 
GCCATGATTACTTCAGCCGT 


(CT) 8 


136 


56 


I 


KF2 18837 


Sebal 3 


F 
R 


VIC-CAAAGCACTAGCATTGCATGA 
GGAAAGGTGCGTTGTTGATT 


(GA) 8 


317 


56 


I 


KF218838 


Sebal7 


F 
R 


6FAM-AGTGGTACTCCGCCAACATC 
TCGTGAATTTGTAGCTTGCG 


(AC) 7 


256 


56 


I 


KF2 18839 


Seba21 


F 
R 


VIC-TTCTTGATGTGGCTCAGCAG 
GCAAACAGAGTTAGCATTGCAT 


(TG) 7 


259 


56 


I 


KF218840 


Seba05119 


F 
R 


NED-CGCCATTTCCAGACAACTCC 
GAGGACGGAGTACGAGAACC 


(CGC) 7 


227 


56 


II 


KF2 18845 


Sebal 1349 


F 
R 


VIC-ATGGGAGCAGGGGTTTACTG 
TGAGCAGGTGGCCTATCATC 


(TGC) 8 


231 


56 


II 


KF218843 


Sebal2491 


F 
R 


ATCGCTTCAGAGCATTGTGG 
AGCCTTGAAAAGGAGATGCC 


(TTG) 7 


248 


56 


NI 


KF2 18847 


Sebal9849 


F 
R 


6FAM-TGTCGCAAAAATGATCTACGTTATG 
CCGAAGAATGGATCGGGAAG 


(ATA) 12 


174 


56 


II 


KF218841 


Seba28414 


F 

R 


NED-GGCTTTCGAACGGAATGGAG 
ACTCCTCCTCAACCCAAACC 


(TTG) 8 


136 


56 


II 


KF218844 


Seba09440 


F 

R 


CCTTACCCTTGTCCCTTCCC 
AGAGACTAATATACTCAGCTTGCG 


(GT) 12 


232 


56 


NI 


KF2 18848 


Sebal 1773 


F 

R 


GAAAGTAGCGGCGGTGTTG 
CTTCAGTCTCAGCTTTGCCC 


(CGG) 7 


247 


56 


NI 


KF218849 


Seba28206 


F 

R 


PET-ACAACATGCAACTTGACCATC 
CTCGTTGCTGGTGATTGAGG 


(CT) 12 


176 


56 


II 


KF218846 


Seba31211 


F 

R 


6FAM-TCCACATCCCATTGGAGCAG 
CTGTGTTTTTAGCATGTGGGAG 


(CA) 15 


250 


56 


II 


KF218842 


Seba31409 


F 

R 


GCACAGTTGTAGGGCATTGAG 
TCCGGTGAATCCTGCTAGAC 


(CA) 15 


241 


56 


NI 


KF218850 


Note: NI 


= not included in multiplex groups; T a = annealing 


temperature. 











a For each forward primer, the fluorescent label is indicated at the 5' end. 
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Table 2. Results of initial primer screening of 12 polymorphic micro satellite loci in three populations of Sebaea aurea: 



Helderberg Nature Reserve (N= 38) Somerset West (N = 25) Paradyskloof (N = 30) 



Locus 


A 


H 0 


H e 




A 




H e 




A 




H e 




oeba(J4 


10 


0.444 


0.874 




1 


0.200 


0.858 


U. / /o^^ 


13 


0.360 


0.917 




SebalO 


5 


0.658 


0.693 


0.052 ns 


5 


0.720 


0.764 


0.059 ns 


6 


0.633 


0.688 


0.080 ns 


Seball 


11 


0.528 


0.847 


0.380*** 


11 


0.571 


0.893 


0.366** 


9 


0.448 


0.844 


0.473*** 


Sebal3 


16 


0.868 


0.775 


-0.122 ns 


9 


0.480 


0.642 


0.256** 


18 


0.833 


0.829 


-0.005 ns 


Sebal7 


6 


0.067 


0.756 


0.915*** 


7 


0.182 


0.806 


0.778*** 


10 


0.296 


0.855 


0.658*** 


Seba21 


9 


0.290 


0.822 


0.650*** 


8 


0.160 


0.777 


0.797*** 


10 


0.621 


0.859 


0.281 ns 


Seba05119 


7 


0.278 


0.824 


0.669*** 


8 


0.286 


0.849 


0.669*** 


8 


0.280 


0.813 


0.660*** 


Seball 349 


6 


0.800 


0.736 


-0.090 ns 


4 


0.583 


0.664 


0.124 ns 


5 


0.633 


0.691 


0.085 ns 


Sebal9849 


2 


0.050 


0.050 


0.000 


1 


0.000 


0.000 




3 


0.069 


0.068 


-0.009 ns 


Seba28414 


5 


0.850 


0.739 


-0.156 ns 


8 


0.792 


0.823 


0.039 ns 


6 


0.500 


0.774 


0.358** 


Seba28206 


4 


0.000 


0.686 


1.000*** 


6 


0.227 


0.610 


0.633*** 


19 


0.571 


0.945 


0.401*** 


Seba31211 


5 


0.750 


0.745 


-0.007 ns 


5 


0.522 


0.693 


0.251 ns 


7 


0.536 


0.684 


0.220* 


Average 


7.17 


0.465 


0.712 


0.353*** 


6.58 


0.394 


0.698 


0.443*** 


9.50 


0.482 


0.747 


0.360*** 



Note: A = number of alleles per locus; F IS = inbreeding coefficient; H e = unbiased expected heterozygosity; H Q = observed heterozygosity; N = number 
of individuals. 

a See Appendix 1 for population locality information. 

b Deviations from HWE were statistically significant at *P < 0.05, **P < 0.01, ***p < 0.001; ns = not significant. 



values. The analysis of progeny arrays would be required to definitely rule out 
or confirm the presence of null alleles. 

The two multiplex reactions were also tested on three to six individuals 
of seven different species of Sebaea (Appendix 1) to assess for potential cross- 
amplification of SSR loci. Cross-amplification was only successful in S. am- 
bigua Cham, and in S. minutiflora Schinz. All 12 micro satellite loci amplified 
in these two species and were polymorphic. Both of these species belong to 
the S. aurea clade, which exclusively contains the tetramerous Sebaea (Kissling 
et al., unpublished data). No successful amplification was obtained in any of the 
five remaining species for the 12 assayed loci. 



CONCLUSIONS 

Twelve microsatellite loci were characterized to investigate the 
reproductive biology and population structure in S. aurea. These 
markers will serve to estimate outcrossing rates in progeny arrays 
and therefore help to understand the reproductive function of 
secondary stigmas in S. aurea. Cross-species transferability ex- 
periments to S. ambigua and S. minutiflora expand the usefulness 
of this set of SSR loci to other species of Sebaea. 
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Appendix 1. Species, collection locality, geographical coordinates, and herbarium voucher of species of Sebaea analyzed in this study. All specimens are 
deposited at the Universite de Neuchatel herbarium (NEU), Neuchatel, Switzerland. 



opecies 




vJcUglapillLal LUUIUlIlalcS 




vuuLiiei iiu. 


N 


S. ambigua Cham. 


South Africa: Western Cape, Pringle Bay 


34°22.627 , S, 18°49.814'E 


40 


Kissling and 
Zeltner 45/2005 


4 


S. aurea (L. f.) Roem. & 


South Africa: Western Cape, Helderberg Nature Reserve 


34°03'43"S, 18°52 , 24 ,, E 


139 


Kissling 23/2010 


38 


Schult. 














South Africa: Western Cape, Somerset West 


34°02'54"S, 18°50'00"E 


212 


Kissling 33/2010 


25 




South Africa: Western Cape, Paradyskloof 


33°58 , 06 ,, S, 018°52 , 43 ,, E 


233 


Kissling 22/2010 


30 


S. exacoides (L.) Schinz 


South Africa: Western Cape, Franschoek Pass along R45 


33°55.041'S, 19°09.518'E 


665 


Kissling and 
Wuerfel 6/2008 


4 


S. marlothii Gilg 


South Africa: Eastern Cape, Naude's Nek 


SO^.U^S, 28°08.394 , E 


2513 


Kissling and 
Zeltner 38/2005 


4 


S. micrantha Schinz 


South Africa: Western Cape, Helderberg Nature Reserve 


34°03.647 , S, 18°52.258 , E 


156 


Kissling 15/2010 


4 


S. minutiflora Schinz 


South Africa: Western Cape, Gordon's Bay, 

unconstructed plot in front of BP station, along R44 


34°08.513 , S,018°51.112 , E 


3 


Kissling 14/2010 


6 


S. spathulata Steud. 


South Africa: Western Cape, Pringle Bay 


34°22.627 , S,018°49.814 , E 


4 


Kissling and 
Zeltner 46/2005 


4 


S. thomasii Schinz 


Lesotho: road to Sehlabathebe lodge 


29°52.692 , S, 029°05.342 , E 


2599 


Kissling and 
Zeltner 26/2005 


3 




Lesotho: ca. 1 km from Schelabathebe lodge 


29°52.058 , S, 029°06.966 , E 


2461 


Kissling and 
Zeltner 29/2005 


3 



Note: N= number of samples analyzed. 
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